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Abstract 

This paper establishes the applicability of density functional theory methods to quantum com- 
puting systems. We show that ground-state and time-dependent density functional theory can 
be applied to quantum computing systems by proving the Hohenberg-Kohn and Runge-Gross 
theorems for a fermionic representation of an iV qubit system. As a first demonstration of this 
approach, time-dependent density functional theory is used to determine the minimum energy gap 
A (TV) arising when the quantum adiabatic evolution algorithm is used to solve instances of the 
NP-Complete problem MAXCUT. It is known that the computational efficiency of this algorithm 
is largely determined by the large- iV scaling behavior of A(N), and so determining this behavior is 
of fundamental significance. As density functional theory has been used to study quantum systems 
with N ~ 10 3 interacting degrees of freedom, the approach introduced in this paper raises the 
realistic prospect of evaluating the gap A(iV) for systems with iV ~ 10 3 qubits. Although the 
calculation of A(iV) serves to illustrate how density functional theory methods can be applied to 
problems in quantum computing, the approach has a much broader range, and shows promise as 
a means for determining the properties of very large quantum computing systems. 



PACS numbers: 31.15.ee,03.67.Ac 



I. INTRODUCTION 



The inability of a classical computer to efficiently simulate the dynamics of a quantum 
system is well-known. The problem is that the dimension of the Hilbert space grows expo- 
nentially with the number of degrees of freedom of the quantum system, and this in turn 
causes an exponential growth in the amount of memory and CPU-time required to carry out 
the simulation. This inefficiency is a major stumbling block for numerical studies aiming 
to determine the asymptotic performance of quantum algorithms. For example, numerical 
simulation of the dynamics of the quantum adiabatic evolution (QAE) algorithm applied to 
the NP-Complete problem Exact Cover 3 has been limited to systems containing N < 20 
qubits 1,2 . Because the algorithm dynamics must be adiabatic, its runtime T must satisfy 
the inequality 



where 



M = max 

0<s<l 



{E 1 {s)\ d ^\E {s)) 



A = mm [Ek{8)-Eo{8)U (2) 

0<s<l 

and here t is time; s = t/T is dimensionless time; H(s) is the time-dependent Hamiltonian 
that drives the dynamics of the QAE algorithm; and {Ei(s), \Ei(s)) : i = 0, . . . ,2 N — 1} 
are the eigenvalues and eigenstates of H(s). In the usual formulation^^ of QAE, dH(s)/ds 
is an s-independent matrix whose largest eigenvalue bounds M. Typically, this eigenvalue 
scales polynomially with N. Thus, if the minimum gap A(iV) separating the ground- and 
first-excited states scales polynomially (exponentially) with N, so will the algorithm runtime 
T(N). An efficient (inefficient) algorithm 6 for a computational problem is one that solves 
all instances of the problem with polynomial (exponential) T(N). We see then that the 
computational efficiency of the QAE algorithm is largely determined by the scaling behavior 
of the minimum gap A(iV). Attempts to evaluate A(iV) using exact diagonalization 7 have 
been limited to iV < 20 qubits. Recently, however, the minimum gap A(iV) for QAE applied 
to Exact Cover 3 has been determined for iV < 128 qubits using quantum Monte Carlo 
methods^. This represents a substantial technical advance, and has stirred great interest in 
finding other computational approaches that might allow quantum algorithm performance 



to be determined for still larger qubit systems. 

Quantum computation is not the only research area struggling with the difficulties of 
simulating quantum systems^. Condensed-matter physicists and quantum chemists have 
been working under the shadow of this problem for decades. A number of computational 
approaches have been developed which, together with increasingly more powerful computers, 
have allowed much progress to be made, in spite of the ultimately unavoidable difficulties 
involved. Among these approaches, density functional theory (DFT) has proven to be one 
of the most successful 1 ^ 1 -* 1 ^. DFT is a theory of interacting fermion systems. It provides an 
exact treatment of all many-body effects through the exchange-correlation energy functional. 
It can also handle the coupling of such fermion systems to both static and time- varying elec- 
tric and magnetic fields. Ground-state density functional theory (GS-DFT) has been used 
to determine a wide range of ground-state properties of atomic, molecular, and solid state 
systems 1 ^ 4 -; while time-dependent density functional theory (TD-DFT) has been used to 
determine excited-state properties, as well as the linear and non-linear response of interact- 
ing many-electron systems to electromagnetic fields 1 ^ 1 ^. For our purposes, it is especially 
significant that DFT has been successfully applied to quantum systems containing N ~ 10 3 
interacting degrees of freedom^* 1 ^ 1 ^. 

In this paper we establish the applicability of DFT methods to quantum computing 
systems. By establishing this link, we shall see that a powerful tool becomes available for 
determining the properties of very large quantum computing systems. Although our analysis 
can be extended to the case of N qudits (o?-level systems) residing on a .D-dimensional lattice, 
we restrict the presentation to N qubits residing on a 2D lattice since this corresponds to 
the experimentally interesting cases of qubits placed in a 2D ion trap^, or restricted to a 
planar superconducting qubit circuit^. 

The outline of this paper is as follows. We begin in Section [Til by showing how an iV qubit 
system can be transformed into a system of N lattice fermions, and then in Section [TTTl illus- 
trate this transformation by using it to re- write the dynamics of the QAE algorithm applied 
to the NP-Complete problem MAXCUT— . For the resulting interacting fermion system, 
Section IIVI establishes the Hohenberg-Kohn- and Runge-Grosa 1 ^ theorems, and sets up the 
auxiliary Kohn-Sham system of non-interacting fermions- 4 . The results of Section ITVl provide 
the justification for applying GS- and TD-DFT to quantum computing systems. The proofs 
given in Section [IV] are adaptations of well-established proofs used for interacting electron 



systems, and so their validity should not be in doubt. Section |V] works out the linear re- 
sponse of the system of interacting fermionized qubits using TD-DFT, and as an application, 
shows how this response can be used to determine the minimum energy gap A(iV) for the 
MAXCUT dynamics. Here we begin to see the value of the newly established link between 
DFT and quantum computing. Calculation of A (AT) boils down to a calculation of exci- 
tation energies, and the reliable calculation of excitation energies for very large interacting 
electron systems was one of the first triumphs of TD-DFT. A straightforward adaptation of 
standard TD-DFT arguments then determines A (A). In light of earlier remarks, the link 
established in this paper between DFT and quantum computing raises the realistic prospect 
of evaluating the minimum gap A(A) for A ~ 10 3 qubits, and thus of studying the perfor- 
mance of the QAE algorithm for much larger qubit systems than is currently possible using 
other approaches. Although we focus on the calculation of the minimum gap in this paper, 
it is clear that the application of DFT to quantum computing systems has a much broader 
range, and shows genuine promise as a means for determining the properties of very large 
quantum computing systems. Finally, the paper closes in Section [VI] with a discussion of 
future work. 

II. QUBIT FERMION TRANSFORMATION 

Consider A qubits residing on an A-site 2D lattice with basis vectors (k = 1, 2), and 
sites specified by the position vector r. Let <x(r) denote the Pauli matrices associated with 
the qubit at r. We now show how the qubits can be converted into lattice fermions via the 
2D Jordan-Wigner (JW) transformation 2 -. Note that the following analysis can be extended 
to A qudits (viz. A d- level systems) on a .D- dimensional lattice using the generalized JW 
transformation 24 that fermionizes a spin s system (d = 2s + 1) in D spatial dimensions. 

For a 2D system of qubits, the JW transformation is: 

<r+(r) = 2alQ r 
<7"(r) = 2Q\a r 
cr z (r) = 2n r — 1. 

Here: c ± (r) = er^r) ± ia y (r); a\ (a r ) creates (annihilates) a lattice fermion at r; n r 



(3) 
= a\.a r 



is the fermion number operator at r; and 

Q r = exp [-i<j) r ] ; r = J $ ( r ' r ') (4) 

In Eq. 03]), <3>(r,r') is the angle made by (r — r') with respect to some reference direction, say 
ei. Thus: (i) $(r,r') changes by 2tt when (r — r') traces out a closed loop around r'; and 
(ii) by convention, $(r, r) = 0. The requirement that the Pauli operators cr(r) commute at 
different lattice sites forces 9 to satisfy 

1 

2m + 1 (m = 0,±1,±2, ■ ••) 



2n9 
in Eq. ®. 

As shown in Ref. 23 



the lattice fermions are spinless, and minimally-coupled to a gauge 
field A k (r) = A k (f> r = 4> r +e k — 4>r- The action for the gauge-field A^{j) is given by the 
Chern-Simons term^ 

J r 

Maxwell's equations for this system take the form 

j? = e^ x F uX (r), (5) 

where is the fermion current, F u \(r) is the gauge field tensor, e fluX is the totally anti- 
symmetric Levi-Civita tensor, and = 0,1,2. From Eq. (jSJ), the fermion current j Tjfl 
has components 



Jr,o = n r , 

jr,k 



^^{A fe G r , y }{9 t% } (k = 1, 2), (6) 
y 

where G r ^ y is the Green's function for the lattice Laplacian 

y~] A fc A fc G riy = -2n5 rt y. 

k=l,2 

Fermion current conservation, d^j^ = 0, follows immediately from Maxwell's equations. 



III. APPLICATION: NP-COMPLETE PROBLEM MAXCUT 



In the problem MAXCUT, one considers an A-node undirected graph Q with nodes 
specified by r. The nodes (edges) are assigned weights w r (w r y), and a binary variable s r 



is associated with each node r. A cut of the graph Q is a partition of the nodes into two 
sets So and «Si. For all nodes belonging to Sq (Si), s r is assigned the value (1). The node 
variables are used to construct a string variable s = s ri ■ ■ ■ s rN , and all possible assignments 
of the N (binary) node variables leads to 2 N possible string assignments for s. It follows 
that there is a one-to-one correspondence between cuts of Q and string assignments for s. 
The MAXCUT problem is to find the cut (viz. string assignment) that maximizes the payoff 
function P(s) given by 



P(s) = S r W r + S r (l — S r r)lV r y. 



MAXCUT is known^ to be NP-Complete, and so it belongs to the set of "hardest problems" 
in the complexity class NP. 



The QAE algorithm was applied to MAXCUT in Ref. |22|, where the dynamics is driven 
by the Hamiltonian 

H(t) = (i-1)h + (I) H P . (7) 



T J V T 
Here T is the algorithm runtime, 

flo = 5>,(r), (8) 

r 

and 

H P = ^ Wr (l-a 2 (r))/2 

r 

+ J2 w ry(l-v z (r)a z ( r '))/2. (9) 

r,r' 

The Hamiltonian Hp is known as the problem Hamiltonian. From Eq. (Q, its eigenstates 
are the simultaneous eigenstates of the {<r 2 (r)}: 

a z (r)\s ri ■ ■ ■ s rN ) = (-l) Sr |s ri •••s TN ). 

By construction^-, each bit string s = Si ■ ■ ■ that maximizes the MAXCUT payoff func- 
tion labels a ground-state |si • • • sjv) of Hp. The QAE algorithm places the qubit system in 
the ground-state of the initial Hamiltonian Hq, and for runtime T sufficiently large, H(t) 
evolves the quantum state adiabatically so that at time T, the system is in the ground-state 
of the final Hamiltonian Hp with probability approaching 1. Measurement of the {a z (r)} at 
time T yields, with probability approaching 1, a string s± ■ ■ -s N that solves the MAXCUT 
instance. 



Using Eqs. ([3]) in H(t) gives the fermionized QAE Hamiltonian for MAXCUT: 



H 



r 




(10) 



where 



w r = v r + W r ; 



and a term proportional to the identity has been suppressed. 

IV. DENSITY FUNCTIONAL THEORY 

In this Section we establish the applicability of the Hohenberg-Kohn and Runge-Gross 
theorems to the QAE/MAXCUT problem. These theorems justify the use of, respectively, 
ground-state and time- dependent density functional theory to the MAXCUT dynamics. 
Throughout, the ground-state is assumed to be non-degenerate, as would be appropriate 
for a non- vanishing minimum gap A. The formalism can be extended, however, to cover 
degenerate ground-states^ 1 ^. 

A. Ground-State Density Functional Theory 

We have seen that the QAE algorithm has an adiabatic dynamics that is driven by a 
slowly- varying Hamiltonian H(t). In this subsection we focus on the MAXCUT Hamiltonian 
H(t) at a fixed instant of time t = t*. By fixing the time, we obtain a well-defined static 
Hermitian operator if* = H(t = t*). The aim of this subsection is to show that the 
Hohenberg-Kohn theorem applies to if*. With this theorem in place, GS-DFT can be used 
to study the ground-state properties of if* = H(t = t*) for any specific intermediate time 
< t* < T. We stress that even though the QAE algorithm works with a slowly-varying 
Hamiltonian, the discussion in this subsection is restricted to the static Hermitian operator 

= H(t*) that is the value of H(t) at the time t = t*. 



Our starting point is the energy functional for the instantaneous MAXCUT Hamiltonian 
if* = H(t = t*): 

E[n] = min (11) 

The domain of E[n] is the set of all iV-representable site occupation functions (SOF) n r that 
can be obtained from an iV-fermion wavefunction. The minimization in Eq. ( TTTT) is over all 
for which 

n = n r = (^\n T \t/)), 

and the minimizing state l^minM) is thus a functional of n r . Let \ip 9 ) denote the ground- 
state of ii*; E 9 the ground-state energy; and n g r the ground-state SOF. Inserting Eq. ffTUj) 
into (TTTj) gives 

E[n] = (U/T)J2 v rn r + Q[n], 

r 

where 

Q[n] = min {^\{T U + U u )\i/>), 

and T u and U tt are the first and third terms, respectively, on the RHS of Eq. ( TTU1) at t = t*. 

To establish the Hohenberg-Kohn (HK) theorem for ii* we must show^: (i) E[n 9 ] = E 9 ; 
(ii) for n r ^ n 9 , E[n] > E 9 ; and (iii) the ground-state expectation value of any observable is a 
unique functional of the ground-state SOF n 9 . By the variational principle, (ip\H*\ip) > E 9 , 
with equality when = \ip 9 ). Thus, for n = n 9 , the search in Eq. (TTTT) returns the 
ground-state \ip 9 ) as the state iVWufn 9 ]) that minimizes E[n 9 ]. It follows that 

E[n 9 } = (^ s |F*|-0 9 ) = E 9 . 

This establishes condition (i). For n ^ n 9 , the minimizing state IVWiM) 7^ IV' 9 ^ 9 ])? an d so 
by the variational principle, 

E[n\ = (Vw[n]|#*|VwM> > E 9 . 

This establishes condition (ii). Finally, since the ground-state \ip 9 ) = \ip m i n [n 9 }), it is a 
functional of n 9 , and consequently, so are all ground-state expectation values: 

(6) gs = (r\o\r) 

= (Vw[n s ]|0|Vw[n 9 ]> 
= 0\n 9 }. 



Condition (iii) is thus established, completing the proof of the HK theorem for H* = H(t*). 

To obtain a practical calculational scheme, an auxiliary system of non-interacting Kohn- 
Sham (KS) fermions is introduced^, and it is assumed that the ground-state SOF n 9 can be 
obtained from the ground-state density of the KS fermions moving in an external potential 
v ks . For if, = H{U), the KS Hamiltonian H ks = T[ + V ks is defined to be: 



H ks = e( 1 -|) fe<4 + fa} + £ (|) 



where q r = (Q r ) is the ground-state expectation value of Q r . The effects of Q r are thus 
incorporated into the KS dynamics through the mean-field q r . The KS energy functional 
e ks [n) is: 

e ks [n] = min (^H^) = T/>] + V (% ) v k r s n r . (12) 
To determine the KS external potential v ks , we re- write Eq. ( II ip as 

E[n] = r u [n] + (|) v r n r + UH (13) 

where 

U[n]=Q[n]-Tl[n] 

is the exchange-correlation energy functional. As noted in Section [U it is through the 
exchange- correlation energy functional £, X c[n] that DFT accounts for all many-body effects. 
Since n 9 r minimizes both €k s [n} and E[n], Eqs. (fl2l and (fT3l) are stationary about n = n 9 . 
Taking their functional derivatives with respect to n, evaluating the result at n = n 9 , and 
eliminating 8T[ /8n\ n =ns gives 

v ? = v *+(£) v ^[n 9 ](r), (14) 

for t* 7^ 0. Here v xc [n 9 ](r) is the exchange-correlation potential which is the functional 
derivative of the exchange-correlation energy functional £xe[n 9 ]: 

SU[n 9 ] 



v xc [n 9 }(r) 



dn- 



This sets in place the formulas for a self-consistent calculation of the ground-state properties 
of = H(t#) using GS-DFT. Entanglement— and its links to quantum phase transitions 30 
have been studied using GS-DFT. 



B. Time Dependent Density Functional Theory 



Here we establish the Runge-Gross theorem^ for the instantaneous MAXCUT dynamics. 
Thus we focus on the instantaneous Hamiltonian = H(t*) for a fixed £* (0 < t* < T). 
Now, however, we suppose that the external potential v r in H(t*) begins to vary at a 
moment we call t — 0. For t < 0, v r (t) = v r , and the fermions are in the ground-state 
\ipo) of H(t*). The Runge-Gross theorem states that the SOFs n r (t) and n' T (t) evolving 
from a common initial state \ip{0)) = \ipo) under the influence of the respective potentials 
V r (t) and Vr(t) (both Taylor-series expandable about t = 0) will be different provided that 
\V r (t) - Vl(t)\ + C(t). For us: 

V r (t) = (|) (l - |) V v {t) 

vm = (|) (i - 1) v'm 

and 

oo 

V T (t) = ^a fc (r)t fc /fc!; 

fc=0 

oo 

k=0 

Let Cfc(r) = afc(r) — a' fc (r). The condition that [V^-(t) — V^'(t)] 7^ C(t) means a smallest integer 
.fT exists such that Cfc(r) is a non-trivial function of r for all k > K, while for k < K, it is 
a constant Ck which can be set to zero without loss of generality. 

Recall (Eq. (jBJ)) that the conserved fermion current has components 

3r,o(t) = n T (t), 

kk(t) = (7^) E(A fc G r>y )^n y (t), 
' y 

with k = 1,2. Defining j T ,k{fy = {^0 |jr,fc(^) IV'o) ? it follows that 

dt{j r ,k(t) = (Mkk(t),H(t) - #'(t)]iVo). (15) 

Here j r ,k(t) (j' r k (t)) and H(t) (H'(t)) are the expected fermion current and the Hamiltonian, 
respectively, when the external potential is v r (t) (v' r (t)). The Hamiltonians H(t) and H'(t) 
differ only in the external potential. Defining 

Sjr,k(t) =j r , k (t) -j' T , k (t), 



and 

5v y (t) = v y (t)-v;(t), 

evaluation of the commutator in Eq. (115)) eventually gives 



d t {Sjr, k (t)} = - £(A fc G r , y ) SV y (t) M y (t), (16) 

y 

where 

M y (t) = (ij \(a y Q y + Q y a y )\ij ). 

With K defined as above, taking K time-derivatives of Eq. (fl6l) and evaluating the result 
at t — gives 



where we have used that 



qK+1 

l(#r,fc(t)) 



(i) E( A ^,y)-M y (0)C x (r), (17) 
' y 



^ (<%(*)) | t=0 = C k (y) = 



for k < K. It is important to note that A4 y (0) ^ 0. This follows since 

[H(Q,n r (Q] ^0 

for ^ T, and so the eigenstates of (specifically, its ground-state l^o)) cannot be 

fermion number eigenstates. This ensures that the ground-state expectation value 

M y (0) = (Mi^Qy + Q y a y )\4j ) ± 
for t*^T. It follows from the continuity equation for the fermion current that 

k=l,2 

Taking K time-derivatives of this equation, evaluating the result at t — 0, and using Eq. (TIT]) 
gives 

QK+2 

(n r (t) - <(t)) | i=0 = -C K (v)M r (0) + 0, (18) 

where we have used the equation of motion for G ry . Equation (Tl8|) indicates that n r (t) 
cannot equal n' r (t) since it insures that they will be different at t = + , and so cannot be 
the same function. This proves the Runge-Gross theorem for the instantaneous MAXCUT 
dynamics. 



We have just seen that when potentials V r (t) and V^(t) differ by a time-dependent function 
C(t), they give rise to the same SOF n r (t). However, the wavefunctions produced by these 
potentials from the same initial state will differ by a time-dependent phase factor. For our 
purposes, it is important to note that this extra phase factor cancels out when calculating 
the expectation value of an operator. In particular, it will cancel out when calculating the 
instantaneous energy eigenvalues E n (t) = (E n (t)\H(t)\E n (t)). As a result, this phase factor 
will not affect our calculation of the minimum energy gap in Section [V] Having said that, 
it is worth noting that this subtlety is not expected to cause difficulties in practice since the 
probe potential V r (t) is assumed to be under the direct control of the experimenter, and so 
the precise form of V r (t) is known. When an experimentalist says a sinusoidal probe potential 
has been applied, this means V r (t) = V r smut; it does not mean V T (t) = V r smujt + C(t). 
Thus in a well-designed experiment C(t) = 0. 

The Kohn-Sham (KS) system of non-interacting fermions can also be introduced in TD- 
DFT 15 . We must still assume that the interacting SOF n r (t) can be obtained from the SOF 
of the non- interacting KS fermions moving in the external potential v^ s (t). The potentials 
v^ s {t) and v r {t) are related via (£* 7^ 0) 



though Eq. ffT9l) is to be thought of as defining the time-dependent exchange-correlation 
potential v xc [n(t)](r) . 

V. MINIMUM GAP 

A problem of longstanding treachery in GS-DFT is the calculation of the excitation 
energies of a fermion system. TD-DFT was able to find these energies by determining the 
system's frequency-dependent linear response, and relating the excitation energies to poles 
appearing in that response. The arguments used^- are quite general, and can be easily 
adapted to determine the energy gap for the instantaneous MAXCUT dynamics. 

Previously, we considered an external potential that becomes time-varying for t > 0. 
Our interest is in the interacting fermion linear response, and so we assume that the total 
potential has the form 




(19) 



v r + vl(t), 



with vl(t) a suitably small time- varying perturbation. The probe potential vl(t) generates 
a first-order response n\{t) in the SOF: 



n t ° t (t) = n 9 T + n 1 T {t). 

The susceptibility \ r y (t — t') connects the first-order probe potential to the SOF response. 
The total potential t>* ot (t) is related to the KS potential v^ s (t) through Eq. (|T9l) . and by 
assumption, the SOF for both the interacting and KS fermions is the same. This allows the 
time- Fourier transform of the SOF response n\(uj) to be determined from the time- Fourier 
transforms of the KS susceptibility x^ s r ,(a;), the exchange-correlation kernel f xc [n 9 ](r, r'; u), 
and the probe potential vl(u): 

y r' 

= £x£M*£M- (20) 

r' 

The KS susceptibility^ depends on the KS static unperturbed orbitals and the corre- 
sponding energy eigenvalues Ej and orbital occupation numbers ff 

XrA") -l^ih- fi) — {i k) + ir] • ( 21 ) 

j,k J 

The exchange-correlation kernel f xc [n 9 } incorporates all many-body effects into the linear 
response dynamics, and is related to the exchange-correlation potential ^[n 9 ] through a 
functional derivative: 

Jxc[ J " 5n9 • 
In general, the interacting fermion excitation energies 

differ from the KS excitation energies 

U jk — Sj — £k- 

The RHS of Eq. (J20l) remains finite as uj — > Qjk, while the first-order SOF response riy(u) 
has a pole at each fl^. Thus the operator on the LHS acting on ny(uj) cannot be invertible. 
Otherwise, its inverse could be applied to both sides of Eq. (1201 with the result that the RHS 



would remain finite as uj — > Qj k , while the LHS would diverge. To avoid this inconsistency. 



the operator must have a zero eigenvalue as u — > flj k . Following Ref. 
following eigenvalue problem: 



31 



one is led to the 



where, writing 



and 



E + = (22) 



Otk'j' — fk' — fj'i 

9* = h(*)M*)> 



we have 



M mj (u) = a k , r J2 ®r> {f xc [n 9 }(v',y';uj)} $ k / 



uj 



It can be shown that X(Qj k ) = 1. 

At this point in the argument, it proves necessary to introduce some form of approxi- 
mation to proceed further. In the single-pole approximation^ the KS poles are assumed to 
be well-separated so that we can focus on a particular KS excitation energy ujj k = uj*. The 
eigenvectors £&y(o;) and the matrix operator M k j- k iji{uj) are finite at uj*, while the eigenvalue 
\(uj) must have a pole there to match the pole on the LHS of Eq. (J22J): 

A( U )= 7 ^L + 0(1). 

[uj — UJ*) 

Let uj* be d-fold degenerate: uj^jx, • ■ ■ ,uJ kd j d = uj*. Matching singularities in Eq. ( 1221) gives 

d 

E M ^;^K) e Kj > = A n (uj*) C kai (uj*), (23) 
i=i 

where i,n — 1, . . . , d. For our purposes, the eigenvalues A n {uj*) are of primary interest and 
are found from Eq. (1231) . From each A n (uj*), we find 

A» = £W . 

(^ - W*) 

Since A n (f^fc) = 1, it follows that the sum of \ n (flj k ) and its complex conjugate is 2. Plugging 
into this sum the singular expressions for \"(Qj k ) and that of its complex conjugate, and 
solving for fi" fe gives 

Q] k = uj* + Re[A n (u*)]. 



Interactions will thus generally split the u;*-degeneracy. Now let 



5E = imnRe[A n (u*)] 

n 

and 

VLjk = minf2" fe . 

n J 

Our expression for f%, then gives 

In the context of the QAE algorithm, our interest is the energy gap 

A(Q =E 1 (Q-E (U) 

separating the instantaneous ground and first-excited states. In this case, our expression for 
Q jk gives 

A(t„) = [e x (Q - e (Q) + 5E(Q. (24) 

To obtain the minimum gap A for QAE numerically, one picks a sufficiently large number 
of U G (0,T); solves for A(t*) using the KS system associated with i?(t*) to evaluate the 
RHS of Eq. f l24l) ; then uses the minimum of the resulting set of A(t#) to upper bound A. 
Because the KS dynamics is non-interacting, it has been possible to treat KS systems with 
iV ~ 10 3 KS fermions^ 1 ^ 1 ^. This would allow the evaluation of the minimum gap A(iV) for 
the QAE algorithm for iV ~ 10 3 . 



VI. DISCUSSION 



As with all KS calculations, the minimum gap calculation requires an approximation 
for the exchange-correlation energy functional £ X c[n]- Note that, because the qubits in a 
quantum register must be located at fixed positions for the register to function properly, the 
associated JW fermions are distinguishable since they are each pinned to a specific lattice 
site. Consequently, anti-symmetrization of the fermion wavefunction is not required, with 
the result that the exchange energy vanishes in the MAXCUT dynamics. The exchange- 
correlation energy functional £, X c[n] is then determined solely by the correlation energy which 
can be calculated using the methods of Ref. 0. Parametrization of these results yields 



analytical expressions for the correlation energy per particle which, upon differentiation, give 
v xc [n] and f xc [n}. Replacing n — > n r in £ xc [n] gives the local density approximation (LDA) for 
GS-DFT; while n — > n r (t) gives the adiabatic local density approximation (ALDA) for TD- 
DFT. These simple approximations have proven to be remarkably successful, and provide a 
good starting point for the minimum gap calculation. Self-interaction corrections to ^ xc [n] 
are not necessary since the two-fermion interaction [see Eq. ffTUj) ] has no self-interaction 
terms. Finally, because the fermions are pinned, it will be necessary to test the gap for 
sensitivity to derivative discontinuities^ 3 - in £ c [n]. 
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